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The formalism for extracting the CKM parameter Vcd from vN production of charm 
is discussed in some detail. The various model assumptions needed are clearly pointed 
out. A direct determination from neutrino induced dimuon production requires vN charm 
■ production data, i/N charm hadronization data, and the semi-muonic branching ratios for 

' charmed hadrons. Hadronization data from FNAL E531 is re-analyzed to take advantage 

of better-determined properties of the charmed hadrons. A small bias in the original 
published result is removed. Neutrino induced charm fragmentation is compared to e'^e~ 
r-^\ charm fragmentation functions; the data are consistent with a common hadronization 

, scheme. An updated value of the mean semi-muonic branching ratio for charmed hadrons 

' produced in vN scattering for > 30 GeV is obtained. This value is used to determine 

Vcd and its associated uncertainties. Prospects for improving the Vcd measurement to 
test the unitarity limit of the CKM matrix are described. 
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I. INTRODUCTION 



' The CKM parameters \Vcd\ and \Vcs\ can presently only be measured via the neutrino production of 

00 ! charm at high energies. This paper summarizes the current state of knowledge of these CKM parameters 

' and estimates the possible precision that will be achieved in future experiments. Much of the material 

I has appeared previously in unpublished form Q , and a condensed summary of this article will appear in 

. a forthcoming review g]. This document is intended to update the previous result and to provide more 

^ ' details for specialists. The timing is propitious as NuTeV, Nomad, and Chorus should shortly produce 

new experimental results. 

^ • In the standard three generation CKM matrix, unitarity and the precise determinations of \Vud\ and 

^ \ \Vus \ tightly constrain \ Vcd\ and \Vcs\- This can be easily appreciated in the Wolfenstein parameterization 

' Q of the CKM matrix, in which the mixing between different generations is parameterized as sin0i2 = A, 

sin023 —A\^, and sin^iae"* = A\'^{p + i-q). In this scheme 

X 



cd 



Vu. 



= 1 + A^X^p + (2.4 X lQ-^)A\p - i), (1) 



and 



^A^A'' ~ 1 - (1.2 X 10-3)yl2. (2) 
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Since A and p are known to be of order one from measurements of Vcb and Vub at CLEO and Argus, | Vcrf| 
and \Vcs\ must be within a few parts per thousand of \Vud\ and |y„s|- 

On the other hand, if three generation unitarity is not assumed the coupling of \Vcd\ to \Vus\ is not 
as tight. For example, in a four generation unitary CKM matrix, mixing between the second and fourth 
generation (sin024) could allow \Vus\ — < 0.03 while maintaining -I- \Vus^ = 1- In summary, 

the standard model predicts \ Vcd\ and \ Vcs\ to a level of ±0.1%. Any larger deviation would indicate new 
physics, and deviations as large as 10% are interesting. 

Section 2 introduces the necessary formalism and points out assumptions that need to be made in 
the analysis of neutrino experiments and external input that is necessary. Section 3 summarizes the 
present state of knowledge of \Vcd\ and \Vcs\ and the auxiliary quantities that are needed to extract 
these parameters. Section 4 considers techniques of improving the measurements and possible systematic 
limitations; and estimates the sensitivity new experiments at Fermilab and CERN. 
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II. NEUTRINO PRODUCTION OF CHARM 



A. Leading Order Quark-Psirton Formalism 

At a fixed reference 4- momentum transfer Qq, the isoscalar cross section for the neutrino production 
of charm can be written to lowest order in QCD as 

d<j{vN IJ.-CX) _ G%ME 

dxdy TT 

(1 - 2^) • - ■ A) ■ [^MEy{l - + - M^] X ' (3) 

{\Vcd? Ql) + dvii, Ql) + us{^, Ql) + dsii, QD] + \Ves\^ 2s(e, Ql)} 

and for anti-neutrinos as 

da{9N fi+cX) _ G%ME 
dxdy TT 

(1 - ■ 0(2/ - 2&) • 2^) • © [2MEy{l - x) + - M^] X • (4) 

{\Vcdf [u{^,Ql) + d{^,Q^)]+\Vosf2s{^,Ql)} 

where in the above two expressions: 

• Gf is the Fermi constant, M is the nucleon mass and E is the incident neutrino energy. 

• TOc is an effective charm mass parameter, and Mc is the lowest mass charmed hadronic system 

allowed by conservation laws. 

• y is the fraction of neutrino energy transferred to the hadronic system, and ^ is the fraction of the 
proton's momentum carried by (massless) struck quark. In terms of the Bjorken scaling variable 

• 5 ['^v{£,,Qo) + dv{^,Qo)] = v{^,Qo) is the valence quark momentum distribution of an isoscalar 
nucleon, and | [us{^,Qo) + ds{^,Qo)] = q{^,Qo) is the light sea quark distribution. Valence and 
sea quark distributions are practically defined via dv{^, Qo) = d{^, Qq) — d{^, Qq) and ds{^, Qq) = 

Qq), and similarly for u-quarks. 

• ■s(C) Qo) ~ ^(?' Qo) ^^'^ strange quark momentum distribution inside the proton. 

The second line in the two cross section formulas above incorporates threshold effects. The first 

factor (1 — 2ME(, ) ^ kinematic factor reflecting the spin i character of the quarks and leptons. The 

2 

second factor Q{y — 2me^ ) imposes the minimum inelasticity requirement for charm production. The 
third factor 6 



^ — 2Mfe imposes the requirement that the neutrino-struck quark system have sufficient 

invariant mass to form a charm quark. The fourth factor B [2M£'y(l — x) + — Mc^] forces the 
invariant mass of the final state hadronic system to exceed the minimum that is compatible with the 
presence of a charmed hadron and a baryon. Real thresholds will not likely be so sharp; and these factors 
are probably approximations. 

Under a number of assumptions detailed below, the major difference between neutrino and anti-neutrino 
production of charm is the possibility of production off the valence quarks in the neutrino case that is 
absent in the anti- neutrino mode. The sensitivity to \Vcd\ follows if one can isolate the valence quark 
contribution to charm production. This may be accomplished in principle by either measuring the charm 
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cross section at high x, where sea quark distributions are smaU; or by subtracting the anti-neutrino 
cross section from the neutrino cross section. The CKM element Vcs always appears in combination with 
s(^, Q^), which is not yet an independently measured quantity. 

The model assumptions that must be addressed are divided into "partonic" and "hadronic" quantities 
below. 



B. Partonic Level Issues in uN Charm Production 

1. Higher Order QCD Effects 

The most obvious question at the partonic level is that of higher order QCD effects. These have been 
addressed by a number of authors Q . While the form of the cross sections becomes much more complex 
in detail, the essential structure remains the same as far as the CKM matrix elements are concerned , 
i.e., 

da{vN fi-cX) _ G'j,ME 
dxdy TT 
iKrfl' Fv [m„ as(g'), f (C, Ql)] + ] (5) 

\Vcd\^ Fs K, as(Q^), v{^, Ql), q{^, Ql), s{^, Ql), G{^, Ql)] + \ 
iKsl' Fss K, as(O'), f (C, Ql), q{^, Ql),s{^, Ql), G(e, Ql)] J 

and 

dajPN fi+cX) _ GlME ^ 

r \V,d\^Fs [fn,,as{Q^),v{^,Ql),q{^,Ql),s{^,Ql),G{^,Ql)] + 1 ' 
I iKsl' Fss K, as{Q^), v{t Ql), q{t Ql), sit Ql), G{t Ql)] j 

In these expressions, Fy, Fs, Fss, Fs, and Fss sue calculable functionals of the parton distributions 
that depend on the running coupling constant as{Q'^), the charm mass rric, and the parton distribution 
functions at the reference Ql . The sea and strange sea functionals Fs and Fss include contributions 
from the nucleon's giuon momentum distribution G{tQl) via PF-gluon fusion. The giuon distribution 
also affects Fs and Fss, but not Fy, via the normal QCD evolution of the structure functions. Higher 
order QCD processes do not affect the dominance of the valence quark distributions at high x; thus |Vcrf| 
is still measurable from the neutrino cross section alone. As long as Fs — Fs and Fss — Fss, the valence 
contribution to charm production can still be isolated by subtracting the anti-neutrino cross section from 
the neutrino cross section. 

These last assumptions seems quite reasonable, however it is not necessary for s(^, Ql) = s(^, Ql), for 

example, at each value of ^, nor even for rf^'s(^', Ql) = d^sl^, Ql). The only rigorous requirement 
is that /p ^s{t ,Ql) — ^s(C';Qo) since the parton distribution functions used here are defined as 
momentum distributions. A mechanism that could generate s{£,,Ql) ^ s(^, Ql) would be the presence 
of a AK~^ state in the proton's wave function Q . If the neutrino interacts with the proton while it is in 
this state, there will be a manifest asymmetry between the s-quark distribution of the A and the s quark 
distribution of the if ^. It would be difficult to combine the neutrino and anti-neutrino data together if 
s{£,,Ql) 7^ Qq) or dsi^Ql) 7^ d{£^,Ql). All that can be stated is that all measurements up to now 
are consistent with these assumptions. 

The most serious theoretical issues in the treatment of higher order terms in neutrino charm production 
are probably the questions of scale dependence and the treatment of mass thresholds in the evolution 
of the running coupling constant. Whether these issues affect CKM matrix parameters depends on 
experimental design. 
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2.. Partonic pr and Mass Effects 



The parton model of deep inelastic scattering assumes that the struck quarks are massless objects 
moving colinear with the proton. Neither of these assumptions need be true. Quarks and gluons inside 
the nucleon will carry transverse momentum, both due to their confinement in the finite volume of the 
nucleon ("intrinsic" pt) and due to hard gluon radiation and qq pair production. These effects appear 
in the longitudinal structure function of the nucleus 

The longitudinal structure function can be calculated perturbatively for high ; at lower values of the 
momentum transfer, higher twist contributions may become important. From the point of view of charm 
production, the most serious question is the effect of the charm quark mass threshold on i?^. 

The charged current charm production cross section contains a longitudinal piece even at lowest order 
just due to the charm mass. To see this simply, consider production off valence quarks, and note that 

1 _ _ -, _ ml/Q^ 

2MEi ~ l + ml/Q'^ ' ^' 

Then, 

da{vdv — > fi^cX) 



dxdy 



1 1 ?/2 
(l-v)-\ ^ H (v- — 

^ ^' l + ml/Q-^ 2 l + m2/Q2 2' 



The decomposition in y isolates the structure functions, from which it follows, neglecting the ^2 term, 
that 

While it seems reasonable to consider the d-quark to be massless, the strange quark might be expected 
to have a current quark mass of wSOO MeV/c^. The modification of the kinematics due to this effect has 
been treated by Tung et al. jij, who have also considered the effects of the target proton mass. There is 
little effect on the differential cross section due to rris or M . 

Other mass effects could enter in a dependence of the Mq threshold on the struck quark type. In the case 
of neutrino- nucleon scattering, the thresholds are Mc = Ma^, Ma^ -I- Af^, and A/ac ^^k depending on 
whether the struck quark is d- valence, d-sea, or s, respectively. Antineutrino charm production lacks the 
valence channel; and the thresholds are slightly higher for d and s struck quarks than for the corresponding 
neutrino case: Mc = Mm + Md for d and Mc = Mm + M^g for s. 

Both parton mass and pt effects are small at large ^ where the valence quark distributions dominate; 
and corrections to the sea quark distributions are the same for neutrinos and antineutrinos. On the other 
hand, higher twist effects are concentrated at high ^, and are thus of concern in a \ Vcd\ extraction. 



3. Nuclear Effects at the Parton Level 



Most measurements of neutrino charm production are made using nuclear targets. One must correct 
for the neutron excess present in the heavier targets. In addition, there may be more subtle nuclear 
effects. 

A standard assumption made in all deep inelastic scattering analyses is that strong isospin symmetry 
is obeyed, i.e., uP{x,Q'^) — d"-{x,Q'^), and dP{x,Q^) = u'"-{x,Q'^); where the superscripts p and n refer 
to proton and neutron, respectively. In addition, it is always assumed that sP(x,Q'^) = s"(x,Q'^). 
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Recent measurements of the Gottfried sum rule by the NMC collaboration |g] indicate that the stronger 
assumption u{x, Q^) — d{x, Q^) is not true at small x. There are then two corrections due to the neutron 
excess in heavier targets, a valence correction for neutrinos only and a sea correction for neutrinos and 
anti-neutrinos: 

= i • {FY K(e, Ql) - dvi^, Ql)] + Ff [uit Ql) - di^, Ql)] } , (8) 

and 

Other nuclear effects are not as obvious; however, so long as all parton distribution functions used in 
charm production analysis are extracted from data using the same or similar nuclear targets, nuclear 
effects should not affect the results at the parton level. Nuclear partonic effects can be important if 
one must take parton distributions extracted from low atomic number targets and use them to make 
corrections in a heavy target. This happens in the case of the neutron excess correction described above. 
The differences between uv{^, Ql) and Q^) are deduced by comparing scattering from hydrogen to 

that of deuterium Q . The results must then be extrapolated to heavier targets by correcting for the EMC 
effect. The EMC effect is determined from electroproduction measurements of the structure F2{x,Q^), 
and not for the valence quark distributions separately. 



C. Hadronization Effects 



1. Charmed Quark Fragmentation 

Experiments can only measure charmed hadron production, not charmed quark production. Some 
experiments, moreover, do not even detect the charmed hadrons; instead only the Icpton (usually muon) 
from semi-leptonic charm decay is measured. Assuming factorization, the charmed hadron (C) cross 
section can be connected to the charmed quark (c) cross section via fragmentation functions: 



da{vN n^CX) da{iyN cX) 

(10) 



dxdydzdp^ dxdy '^'t)' 

da{9N ^ H+CX) _ da{i?N ~* ^i+cX) 



dxdydzdpTp dxdy 

Here, D'^{z,p^) is the probability distribution for the charmed quark fragmenting into a charmed hadron 
of type h carrying a fraction of the quark's longitudinal momentum z and transverse momentum px with 
respect to the quark direction. The number fh is the mean multiplicity of the hadron h in neutrino 
production of charm. The analogous objects for anti-neutrinos are indicated by the barred quantities. 
Since only one c-quark is produced in a charged current interaction, one can set the normalization 
conditions as 

dz dp'^D^{z,p'^) = dz dplD'^iz.pD^l, (11) 
Jo JO Jo Jo 

and 

E^ = E^ = i- (12) 
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2.. Neutrino vs. Anti-neutrino Hadronization 



It need not be the case that fh = fh or D'^{z,p'^) = D^{z,p^) smce the remnant nucleon will not be 
the same in vN and DN scattering at low energy. The threshold behavior, for example, differs: 

but 

At sufficiently high energies, one expects the charm quark hadronization to become independent of the 
remnant nucleon, and consequently I?J^(z,p^) — > D^{z,p'^) and fh — fh 



3. Quasi-elastic Charm Production 

The cross section for i^n /i^Ac,/i~Sc were first calculated by Shrock and Lee Using SU(4) 
flavor symmetry to relate a{i'n — )■ /i^Ac) to (T(j/n — > , these authors calculated a cross section of 

a{i'n Ac) — 2.4 fb for E^, > 20 GcV. The Fermilab E531 experiment ||l^ observed three events 

consistent with quasi-elastic charm production, from which they obtained a cross section of (O.37I023) 
fb. The overestimate of the theoretical calculation is likely due to the substitution of charmed vector 
meson masses for the p and Ai masses appearing in the vn — > fj,~p cross section. Converting the E531 
data into the upper limit a{m fJ-^^c) < 1.0 fb at 90% C.L., and assuming that the cross section 
is independent of energy, one finds that the fraction of neutrino charm production that appears in the 
quasi-elastic channel is less than 88%, 24%, 13%, 6%, and 2.5% for =10, 20, 30, 50, and 100 GeV, 
respectively. These numbers allow one to judge the extent to which the cross section factorizes into 
independent production and hadronization terms. Factorization is clearly not valid for E^ < 20 GeV. 

Other non-quark-parton model calculations of neutrino charm production exist. Einhorn and Lee |^ 
considered a generalization of vector and axial vector meson dominance (VDM) to the charm sector in 
the same spirit as Shrock and Lee's Ac calculation. The VDM contribution to the neutrino charm cross 
section calculated in this model is very large, due again due the replacement of light vector meson masses 
with charmed vector masses in the relevant form factors. The agreement of the data with the quark 
parton model of charm production at high energies and the failure of the VDM to describe the non- 
charm cross section indicates that VDM is not a dominant contribution at high energies. For E^, < 20 
GeV, this mechanism could well be important. 



4- Neutrino Induced Di-lepton Production 

Emulsion experiments can, in principle, identify all charmed hadrons via their finite decay length. 
In practice, there will be efficiency variations according to the charged multiplicity of the modes; and 
acceptance for low z will be limited, particularly at low vales of E^, or y. 

Other experiments that detect only leptons must de-convolve another layer from the decay process: 



dxdydzdp^dk*d cos 9* dxdy 



E^-^'(^,PT)®s;:T;:(r,cosr), (i3) 

h 

with B'll the semi-leptonic branching fraction for the charmed hadron h, and r^(fc*,cos0*) the joint 
distribution function for the muon momentum in the hadron rest frame k* and the decay angle of the 
muon with respect to the hadron direction in the hadron rest frame. The '(8)' symbol indicates that a 
boost must be performed along the charmed hadron direction to get to the hadron rest frame. To go 
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from a dimuon measurement to a charm production measurement thus requires knowledge of the relative 
production of the various charmed hadrons, their fragmentation functions, their semi-leptonic branching 
fraction, and their decay distributions. 

It is possible that the nuclear environment could modify the charmed quark fragmentation functions 
for production off heavy targets relative to those measured in e^e~ experiments. Recent measurements 
[ pl| of D meson production using hadron beams on nuclear targets are consistent with an A^ '^ behavior 
of the total production rate and of the xf and pT distributions. This is evidence against the presence of 
nuclear effects in either the production or fragmentation of charmed quarks. 

Because of the fragmentation and decay parameters present, neutrino di-lepton experiments can only 
measure the quantity • |T4dP, where the mean semi-muonic branching ratio, defined via 

B,=Y^h-Bl, (14) 

h 

is assumed to be the same for neutrinos and antineutrinos. A calculation of B^^ requires knowledge of both 
the production fractions fh and the semi-muonic branching fractions B^. In addition, the fragmentation 
functions D'^{z,p^) and the decay distributions must be available in order to correct for acceptance and 
smearing. In practice, a mean fragmentation function to all charmed hadrons is used for this task. 



D. Methods for Extracting \Vcd\ and \Vcs\ 

The most complete analyses for the extraction of the CKM parameters in neutrino scattering are those 
of the CCFR collaboration. Distributions of "visible" versions of the kinematic variables E, x, and z 
described above are fitted to a model with four parameters: 

• the charm mass; 

• a, a parameter that describes the difference in shape between s(^, Q^) and q{^, Q^); 

• K, the ratio of the nucleon's momentum carried by s quarks relative to that carried by u and d sea 
quarks, defined as 

..4^, (1.) 

• Bfj^, the mean semi-muonic branching ratio for charmed hadrons produced in neutrino anti-neutrino 
scattering. 

The two analyses fit the same data to the cross section described in equation |l3 using either leading 
order or next-to-leading order [|l3| QCD. To extract values for k and B^, \Vcd and \Vcs\ are set to 
values implied by three generation CKM unitarity. If the unitarity constraints are dropped, then the 
CCFR fit parameters change from B^ ■ \Vcd\'^ and k ^B^ ■ -^p^ |V^s|^- 

The CCFR analyses use the maximum information in each dimuon event. However, the procedure 
couples the desired CKM matrix elements to the properties of the strange sea. Furthermore, the fits to 
absolute rates are sensitive to QCD scale uncertainties. 

A cleaner way of obtaining \Vcd\ from data is to measure the ratio of differences 



This can be recast in a more useful form for experiments: 
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(17) 



where r^t^ = '^^ai^N-'^fi- x)"^ ''mm ~ ^a(t>N^^i+ x) ^ ^'"^ normalized dimuon rates, and r 

is the well-measured D Iv chare;ed current cross section ratio. 
To leading order in QCD, 



3^ \y ,2 K[m,,E,v{£,.Q^)] 



where 



ifK,£;,i;(e,Q2)]^ 

djdy^X ^ ^)9(y - ^)9(^ - ^)Q{2MEy{\ - x) + ~ Mg,]z;(g, Q^) (19) 

/dez;(e,Q2) 

The actual experimental procedure is to measure as a function of neutrino energy. A fit can then be 
performed to the data with T4dP and vHc the free parameters. One still must get the mean branching 
ratio from other data to obtain the CKM matrix element. 

In an emulsion experiment, one could improve the situation by measuring directly the ratio 

" ~ (j{vN n-X) - a{DN fi+X) ' ^ ' 

Again, this can be expressed as a ratio of differences 

= (21) 

with Tc = '^^^^^Z'^^-x) ^^"^ ~ "a{^N^^i+x) • leading order expression for is the same as for 
R~^ with the substitution 1. 

The nice property of R~^ and R^ is their functional dependence only on the valence quark momentum 
distribution. This fact has three virtues: the valence quark distribution is better measured; the valence 
quark contribution to charm production suffers less suppression from the charm mass; and the valence 
quark distribution has a theoretical simpler QCD evolution. Also, since R~^ depends only on ratios of 
cross sections, it should be less susceptible to QCD scale errors. The cleaner systematics is offset by the 
need for high statistics since R~,, is a ratio of differences. 

^ MM 



III. CURRENT STATE OF jVcol AND |Vcs| 

A. Measurements of ■ \ Vcdf 

The most up-to-date extraction of i?^ • \Vcd\^ is from the CCFR collaboration ||l3|. Employing a 
next-to-leading order QCD formalism: 

Bf, ■ \Vcd\^ = (5.34 ± 0.39 ± 0.24l°J5) x lO'^ (CCFR - NLO). 

The first two errors are the experimental statistical and systematic errors, respectively, the systematic 
error being dominated by the uncertainty in modeling the charmed quark fragmentation into charmed 
hadrons. The final error is attributed to the QCD factorization and renormalization scale uncertainties. 
This analysis also indicates that there is a negligible difference in the value of i?^ • jVcdP if the analysis 
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is performed to leading order or next-to-leading order in QCD. Accordingly, one may also use the older 
result from the CDHS collaboration |l4[ : 

■ iKdP = (4.1 ± 0.7t°o,ll) X 10-3 (CDHS - LO). 

The first error is the total experimental error. The second error is the QCD scale error, which is not 
given by the original analysis, but is instead assumed to be the same as in the CCFR measurement. 

Combining the two results, assuming that all of the experimental errors are completely uncorrelated, 
but that the QCD scale error is totally correlated and equal to the CCFR value, yields 

5^ • \Vcd\^ = (5.02l°;^0) X 10-3 (CCFR/CDHS). 

The CKM element Vcs is contained in the CCFR measurement 

• \Vcsf = (2.00 ± 0.10 ± o.oet":?^) x lo-^, 

where the uncertainties represent statistics, experimental systematics, and the QCD scale, respectively. 

B. Charm Production Fractions 

The charm fractions fh have only been measured directly in one experiment, FNAL E531 pc| ]. In 
checking over the E531 result, a bias was detected in the way that they extracted their charmed hadron 
production fractions. Their data is re-fitted with the bias removed. If one believes that fragmentation 
functions are universal, one can check the re-fitted E531 measurements against similar fractions measured 
in e+e~ experiments with similar kinematics. The results of this check are also presented. 

1. Reanalysis Of E531 Data 

In E531, charm could be tagged by the presence of a detached secondary vertex in the emulsion target. 
Some 122 events were tagged in this way, 119 which are neutrino induced and three induced by anti- 
neutrinos. The charge of each charmed hadron was determined by counting the number of prongs coming 
from its decay vertex. Under the assumption that the production of neutral charmed baryons is negligible 
at E531 energies, all neutral candidates are unambiguously identified as D° mesons. (There is one neutral 
particle in the sample that had an identified proton; this event is a candidate for a neutral charmed baryon, 
most likely the a csd state. However, the reconstructed lifetime of the track is very long, indicating 
that this event could be background, or that it could have be kinematically mis-fitted.) Charged charmed 
particles, on the other hand, could be one of three possibilities (again neglecting heavier baryons): a , 
a Dg, or a A J. A good fraction of the Ajt could be identified by the presence of a proton in the emulsion. 
The majority of the remaining charged events could be a or Dg with essentially equal probability. 
The bias in the published E531 result is that they resolved this ambiguity by counting all of the "toss-up" 
events as mesons. They accounted for the possibility of a mistake in this procedure by increasing 
the error in the charmed fractions. Given the state of knowledge of charmed mesons at the time, this 
procedure was not unreasonable. 

2. Re-fittmg E531 Data 

Since the complete data set from the experiment is readily available in the thesis of S. Frederikson 
one can re-do the analysis to remove the bias. The procedure is to try to use all the available 
information about each event to construct a likelihood function. The relevant items are the relative 
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kinematic fit probabilities to the vs. and the proper hfetimes (assuming a particular mass 

hypothesis). Unfortunately, the kinematics offers essentially no separation because of the small 13^, £)+ 
mass difference. The lifetimes, on the other hand, are quite different for the two mesons. Accordingly, a 
likelihood function is constructed for each event using only the decay length information. 

The form of the probability function for each event with a charged charmed hadron is taken to be 

EM;£(7l) 

nn) - ' . (22) 

^ — 

In this expression: 

• The index i takes on values Z)+, Dg, A^. 

• The Ni are the number of produced hadrons of each species type. These are the free parameters in 
the fit. 

• The measured quantities £{n) and qi(n) are the decay length and the momentum for the kinematic 
fit to the hadron i , respectively. 

• The Mi are the masses of the charmed hadrons. From the 1992 PDGR: Mdo = 1.8645±0005 GeV/c^, 
Md+ = 1.8693 ± 0.0005 GeV/c^, = 1.9688 ± 0.0007 GeV/c^, and = 2.2849 ± 0.0006 

GeV/c2. 

• The Ti are the mean lifetimes of the particles. From the 1992 PDG: r£,+ = (10.66 ±0.23) x lO^^^s, 
Too = (4.20 ±0.08) X 10-"s, Tjj+ = (4.50 ±0.28) x lO^^g, and t^h- = (1.91 ±0.14) x lO^^s. Note 
that these are known much more accurately than when E531 ran. 

• The Wi are the acceptance weights for each hypothesis and each event. 



= Qiin) 



di-Mi ^ 



-e "^i-'li ■ t{l) 



(23) 



These weights include the effects of the finite size and resolution of the emulsion, of a special pt 
cut for one prong charm events, and for general detector acceptance. Frederikson's thesis gives the 
weight only for one hypothesis for each event. Weights for the other hypotheses are estimated by 
assuming that the dominant difference in the decay weights is due to the interplay of the lifetime 
and the minimum resolvable decay distance. If the minimum resolvable decay distance is tmim 
then, approximately. 



~ e i^^i , (24) 

so that 

log(wj) Mjq^n ' 



^Most of these measurements have now been improved in the 1996 PDG summary; however the relative change 
between the 1992 and 1996 PDG is very small compared to that between the E531 publication and the 1992 
values used here. 
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independent of imin- Thus, if the assumption is correct, one can scale the acceptance weights using 
the one given weight and the Ufetimes, masses, and momenta. This scahng should take into account 
more complicated effects like the px cut placed on one prong events. An explicit calculation using 
a minimum cut-off of irnin = 15 fim and an analogous maximum cut-off of £max =1.5 cm was also 
tried, with no change in the results. Qi{n) is set to if some information about the event rules 
out a particular hypothesis. This occurs typically if no good kinematic fit existed for one of the 
hypotheses, or if the presence of an identified proton in the final state required the charmed hadron 
to be a baryon. Otherwise Qi{n) ~ 1. 

From the probability functions for each event, an extended log-likelihood function is constructed: 

i = - ^ logP„ - Nobs \0g{NDWD + NdsWDs + ^^Ac^Ac) + {NdWd + NdsWDs + Na^WAc) 
n 

The second and third terms above are the log of the Poisson probability function for observing Nobs 
events given the produced events. The bars over the acceptance weights indicate that these quantities 
are averaged. The Poisson term incorporates the finite statistics of the experiment. The (negative) 
log-likelihood function is then minimized. The results of fits are given in the next section. 



3. Results of Fits 

The fit results are given in Table | in the form of charmed hadron fractions: /do, fo+, /d+i ^''^d /ac- 
The fractions, their errors, and their correlations are obtained from MINUIT. (As a technical aside, the 
Ni defined in the previous section are expressed in terms of the fi and the total number of produced 
charmed hadron events Nq- The fi and Nc are then allowed to vary. The constraint '^2 fi — 1 keeps the 
number of free parameters the same.). The correlations between the fractions are give in Tabic p. 



Energy (GeV) 




+ 




/a+ 


5-20 


0.32 ±0.11 


0.05 ± 0.06 


0.18 ±0.10 


0.44 ±0.12 


20-40 


0.50 ±0.08 


0.10 ±0.08 


0.22 ± 0.08 


0.18 ±0.07 


40-80 


0.64 ± 0.08 


0.22 ±0.09 


0.09 ± 0.08 


0.05 ± 0.04 


> 80 


0.60 ±0.11 


0.30 ±0.11 


0.00 ± 0.06 


0.09 ± 0.08 


> 40 


0.61 ±0.06 


0.27 ±0.03 


0.04 ±0.01 


0.07 ±0.02 


> 30 


0.58 ±0.06 


0.26 ± 0.06 


0.07 ±0.05 


0.07 ±0.04 


> 20 


0.56 ±0.05 


0.20 ± 0.05 


0.11 ±0.04 


0.11 ±0.04 


> 5 


0.53 ±0.05 


0.16 ±0.04 


0.13 ±0.04 


0.17 ±0.04 



TABLE I. E531 Re-fitted Production Fraction Results 



Energy (GeV) 


C(D",D+) 








C(O-^AJ) 


C(i3jA+) 


5-20 


-0.165 


-0.332 


-0.572 


-0.116 


-0.227 


-0.451 


20-40 


-0.166 


-0.280 


-0.229 


-0.810 


±0.034 


-0.332 


40-80 


-0.497 


-0.302 


-0.239 


-0.557 


-0.082 


-0.205 


> 80 


-0.633 


-0.001 


-0.304 


-0.002 


-0.410 


-0.001 


> 40 


-0.619 


-0.166 


-0.288 


-0.382 


-0.250 


-0.150 


> 30 


-0.503 


-0.215 


-0.268 


-0.551 


-0.159 


-0.192 


> 20 


-0.404 


-0.284 


-0.375 


-0.544 


-0.181 


-0.138 


> 5 


-0.344 


-0.329 


-0.433 


-0.466 


-0.184 


-0.208 



TABLE II. Correlations Among Production Fractions 
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4- Theoretical Expectations 



One can make a simple model for what one expects the charm branching fractions to be: 



V+l ' 

= (1 - A) • (1 - • +;f "^ 

/c+ = (1 - A) • 0; 
' /a = A. 

The parameters above have the following meanings: 

• A is the fraction of time the c-quark fragments into a . 

• (f> is the relativity probability the c-quark fragments into a Dg vs. any charmed meson. 

• V represents the relative probability of a c-quark fragmenting into a vector D* meson compared to 
the probability of fragmenting directly into a pseudoscalar D meson. 

• The are D* branching ratios. From the 1992 PDG: B^^ = B{D*" ^ D"X) = 1.00, = 
B{D*° D+X) = 0.00, B++ = B{D*+ D+X) = 0.45 ± 0.04, 5+° = B{D*+ D°X) = 
0.55 T 0.04. 

On the basis of spin counting rules, one expects V = i. Studies of non-charm fragmentation in e+e^ 
experiments indicate that (f) « 0.15 and A w 0.10. The model thus predicts f^o = 0.54, = 0.22, 
fp+ = 0.13, and /a = 0.10. 



5. e+e" DATA 

The fractions /^can be obtained from CLEO data |l^ if one accepts the hypothesis that fragmentation 
functions obtained from e^e~ experiments with ^ — 10.55 GeV can be used to describe quark fragmen- 
tation in neutrino induced interactions with mean final state hadronic mass < W 10 GcV/c^. CLEO 
measured the product of branching ratio times cross section, a ■ B, for charged and neutral D and D* 
mesons, Dg mesons, and A^t baryons for z = > 0.5. They then extrapolated using a tuned Lund 

fragmentation model to all z. Its is possible to convert their a ■ B values into cross sections using up-to- 



date charmed hadron decay rates from the 1992 PDG. Table III summarizes the CLEO measurements of 



(T • B and the branching ratio corrected cross sections integrated over all z and for z > 0.5. 



Mode 


G ■ B(pb) 


cr • B (pb) 


B.R.(%) 


a (nb) 


cr (nb) 




2 > 0.5 


all z 




2 > 0.5 


all 2 




27.0 ± 1.4 


52 ±6 


3.65 ±0.21 


0.74 ±0.06 


1.42 ±0.19 


D+ K-TT + TV+ 


29.3 ±3.0 


47 ±7 


8.0 ±0.8 


0.37 ±0.05 


0.59 ±0.11 


D*+ -^■K+D'\K-n+) 


10.9 ± 1.0 


17 ± 2 


2.00 ±0.18 


0.54 ±0.06 


0.85 ±0.13 


D*+ ~* tv+D'^Ik-tv+tv+tv-) 


23.1 ± 1.0 


33 ±3 


4.1 ±0.3 


0.56 ±0.06 


0.80 ±0.12 


D*° ^'y/pi°D°{K-7v+) 


19.8 ±4.0 


30 ±7 


3.65 ±0.21 


0.54 ±0.11 


0.82 ±0.19 


D+ 4)%+ 


5.8 ±1.0 


7.2 ±2.0 


2.8 ±0.5 


0.21 ±0.05 


0.26 ±0.09 


Aj ^ pK-Tv+ 


8.6 ±1.0 


13.5 ±4.0 


3.2 ±0.7 


0.27 ±0.07 


0.42 ±0.16 



TABLE III. Charmed Hadron Production Measurements By CLEO 
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From the data in Table [11 one can extract th e fra ctions fi as measured in e~^e scattering atyG — 10.55 
GeV. These fractions are summarized in Table IV 



fraction 


z > 0.5 


all z 




0.47 ±0.03 


0.53 ±0.05 


fD + 


0.23 ± 0.03 


0.22 ±0.03 


^Dt 


0.13 ±0.02 


0.10 ±0.06 




0.17 ±0.04 


0.16 ±0.05 



TABLE IV. Charm Hadron Fractions Measured By CLEO 



13 



Using the D* and D data together, CLEO also determined 



= 0.85 ±0.11 ±0.17. 

V + 1 

This value is consistent with the model expectation of 0.75. 

Electron-positron data thus supports the simple model presented above and yields charmed hadron 
fractions in agreement with the re-fit E531 results, but not with the published E531 values. 



C. Updated Charm Branching Ratios 

1. Direct Extraction of 

The charm branching ratios were measured in the early to mid eighties at SLAC. The measurements 
are not precise; one, the £>g semi-muonic rate, is essentially not known at all. One can attempt to 
use a weak form of the spectator model of charm decays to improve the precision, but the correlations 
introduced result in no practical gain. 

The semi-muonic branching ratios for most charmed particles have changed little recently, with one 
important exception. (Actually, only one measurement exists for the semi-muonic rate; the remainder are 
semi-electronic rates that are converted assuming e — universality.) All values below come from the 1996 
PDG except for the D", which also includes a substantially improved number from CLEO [^Sj, yielding 
Bf,{D°) = (6.75 ± 0.30)%. The D+ rate is essentially that of the Mark III collaboration: B^{D+) = 
(17.2±1.9)%. TheDj rate is from a single poor measurement by Mark III ||l|: B^{D+) = (5.0±5.4)%. 
The rate is an ancient measurement from Mark II (at SPEAR) B^{K^ ) = (4.5 ± 1.7)%. 

Using the direct measurements alone and the re-fitted E531 data yields 

B^ = 0.0919 ± 0.0085cF ± 0.0041br , 

where the first error is the contribution due to the charmed hadron species fractions and the second error 
is due to the charmed hadron semi-muonic branching ratios. Charm hadron production fractions now 
dominate this error. 



2. Spectator Model Fits to Branching Fractions 

The I?g and A^t, branching ratios are poorly measured. One can try to improve the situation by 
imposing a weak form of the spectator model which requires the semi-electronic partial width Fe be the 
same for all charmed particles. The semi-muonic branching fraction can then be related to the lifetime 
of a charmed hadron Xn via 

One can use the spectator model to fit for improved values of the branching ratios by minimizing 

2 



2 IB 

B',] - B'^,] {'rt^'^h 



with respect to improved values of the branching ratios B^ and the common semi-electronic width Eg. 
The results of the fit are: ^^(D") = (6.76±0.30)%, ^^(13+) = (17.2±0.8)%, By.{D%) = (7.6±0.4)%, and 
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Bfj,{A'^ ) = (3.4 ±0.3)%. The spectator model fit introduces large correlations among the fit parameters: 
Cdod+ = 0.930, C^oc+ = 0.613, C^oaJ = 0.698, C^+^+ = 0.598, C^+^J = 0-681, and C^+a+ = 0.449. 
The correlations must be included in calculating the error on the fitted mean branching ratio B^. The 
result is = 0.0930 ± 0.0088. 



D. The CKM Parameters 



1. I Vcd I from Neutrino Scattering 
Using the updated value for obtained above: 

\Vcd\ = 0.2S2t°aml (DIRECT). 

The CKM parameter \Vcd\ is thus known to ±9% from direct measurement. Of this, the contribution 
from the uncertainty on is ±6.6%, from the experimental error on B^ ■ \Vcd\'^ 4.2%, and from the 
QCD scale uncertainty t5'3%. The measurement is consistent with the value inferred from the assumed 
unitarity of the CKM matrix of 

\Vcd\ = 0.221 ± 0.003. (UNITARITY), 
but the unitarity prediction is clearly not being tested. 



2. Other Measurements of \Vcd\ 
There are no other direct measurements of I Vcd I- The closest are the measurements of 



Vcd 


2 


fm 


2 


B{D - 


ttIv) 


Vcs 




/^(O) 




B{D - 


^ Klv) 



by Mark III in the neutral D mode and CLEO in the charged D mode |2J]. In this expression, 
B{D — > Triv) and B{D — s- K^v) are the measured Z?£3 semi-leptonic branching fractions; e is a precisely 
known kinematic factor; and /+(0) and f^{Q) are the transition form factors evaluated at zero 
momentum transfer to the system. The Mark III measurement is 



= 0.057+0 0^5 ± 0.005(MarkIII). 





Vcd 


2 


/^(o) 




Vcs 




/^(O) 


The CLEO measurement is 










Vcd 


2 


/-(O) ' 




Vcs 




/f(0) 


Combining the two results yields 






Vcd 


2 


/;(o) 




Vcs 




/f(0) 



0.085 ± 0.027 ± 0.014(CLEO). 



= 0.069 ± 0.020(COMBINED). 



The ratio 



/:(o) 



/f (0) 



is model dependent and ranges from 0.7— 1.4 [2^. This translates into 30% uncertainty 

in \Vcd\- If the theoretical uncertainty in the form factor ratio could be reduced to 10%, then the error on 
\Vcd\ from D meson decay would be ±0.034 or 15%. This is less precise than the neutrino measurement. 
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3. Limit on the Wolfenstein Parameters A and p 



The quantity 



1 can be used, via equation |l|, to set the constraint 

< 173 at 90% C.L. 



A^- 



P- 



This is obviously not a very exciting hmit. There is no evidence for the third generation from quark 
mixing measurements involving only the first two generations. 



4- Limit on \Vcs\ 



If one assumes that the strange sea carries no more momentum than the light quark sea, i.e., k < 1, 
then it follows from the CCFR next-to-leading-order analysis that 

\Vcs\> 0.74 at 90% C.L. 

If one arbitrarily assumes that k, = 0.5 ± 0.5, then 

\Vcs\ = 1.04 ±0.42, 

and 

Vcd 



= 0.23 ±0.09. 



The errors in these latter two quantities is dominated by the uncertainty in k. Note that the ratio 
IVcd/V'csl derived from neutrino and antineutrino scattering is independent of -B^. 
If K can be independently measured to an accuracy of ±25% from the structure function difference 



xF^ — xF^, then the errors on IV^^I and 



Vcs 



would be reduced to 10%. 



IV. FUTURE MEASUREMENTS OF \Vcd\ IN lyN SCATTERING 

Four high energy neutrino experiments are either now running or are approved to run in the next six 
years. While none of the experiments are optimized for the study of neutrino charm production, all have 
the potential to improve the CKM matrix element measurements. The experiments are summarized in 
Table 0. The Nomad |^ and Chorus experiments at CERN are designed to search for j/^ Vr 
oscillations. Nomad features a low mass target with very good tracking and electron identification. This 
experiment should be able to detect charm in both di-lepton modes (/i/x and pe). Their excellent tracking 
may also allow for the identification of charm via the "I?* Dtt trick" . Chorus is a hybrid emulsion 
spectrometer. It's major virtue is its ability to reconstruct charm inclusively via the identification of the 
charm decay vertex. This feature serves to boost statistics, and, more importantly, largely eliminates the 
need to know the production, fragmentation, and decay properties of the charmed hadrons. Fermilab 
E815 uses the E744/770 Lab E neutrino detector. The experiment is optimized for precision studies of 
neutral current interactions. The feature most relevant for charm studies is the new sign-selected neutrino 
beam. This will eliminate the v/P confusion in the dimuon channel and permit a cleaner measurement 
of \Vcs\, assuming that the strange sea is independently known by then. 

The uhimate neutrino charm production experiments are COSMOS (FNAL E803) ^ and TOSCA 
[0 at CERN. Like Chorus, these experiments are designed for a high sensitivity search for i^^ v^- 
oscillations using hybrid emulsion spectrometers. E803 will have a factor of twenty higher statistics than 
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Chorus; and its spectrometer will have three times better resolution. The higher resolution is crucial to 
reduce backgrounds, particularly in one-prong decays of charm. E803 might be able to achieve a resolution 
of ~ 2% on \Vcd\- This is estimated by assuming: a sample of 50,000 reconstructed charm events, which 
reduces the statistical error to ±0.003; a x5 reduction in the experimental systematic errors due to the 
elimination of fragmentation uncertainties and background; a x5 reduc;tion in the QCD scale error via 
the normalization of charm to single muon production that is possible with higher statistics; and a x 10 
reduction in production fraction and branching ratios achieved by the ability to inclusively reconstruct 
charm. The total ±0.004 error on |Vcd| will be comparable to that on \Vus\', and one will thus be able to 
test the unitarity property of the CKM matrix at a level that is sensitive to new physics. 



Experiment 


Target 


Start 


CC Sample 


Charm Sample 


Nomad (CERN) 
Chorus (CERN) 
NuTeV (FNAL) 
COSMOS (FNAL) 


low mass 
emulsion 

iron 

emulsion 


1994 
1994 

1996 
2001 


1 X 10** 
3 X 10^ 

3 X 10" 
8 X 10'^ 


2 X W^ieiJ^ufi, D*) 
2 X 10'' (inclusive) 
2 X 10* (^^) 
4 X 10* (inclusive) 



TABLE V. Future Neutrino Experiments. Event samples are rough estimates. 
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